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^ Summary. The peculiar velocities of galaxies distort the pattern of galaxy clustering 

in redshift space, making the redshift space power spectrum anisotropic. In the linear 
regime of gravitational instability models, the strength of this distortion depends only 
on the ratio fi = f(^l)/b ~ OP'^ jb, where O is the cosmological density parameter 

^) and b is the bias parameter. We derive a linear theory estimator for [i based on the 

harmonic moments of the redshift space power spectrum. Using iV-body simulations, 
we examine the impact of non-linear gravitational clustering on the power spectrum 

Q>^ anisotropy and on our /3-estimator. Non-linear effects can be important out to wave- 

lengths A ~ 50h~^Mpc or larger; in most cases, they lower the quadrupole moment of 
Q_i the power spectrum and thereby depress the estimate of [i below the true value. With 

^ a sufficiently large redshift survey, the scaling of non-linear effects may allow separate 

determinations of O and b. 

We describe a practical technique for measuring the anisotropy of the power spec- 
trum from galaxy redshift surveys, and we test the technique on mock catalogues 
drawn from the iV-body simulations. Preliminary application of our methods to the 
1.2 Jy IRAS galaxy survey yields /3est ~ 0.3 - 0.4 at wavelengths A ~ 30 - 40h-iMpc. 
Non-linear effects remain important at these scales, so this estimate of /3 is probably 
lower than the true value. 



1 Introduction 

A fundamental assumption of the standard cosmological model is that the universe has no pre- 
ferred direction, so that galaxy clustering is, statistically speaking, isotropic. However, we usually 
map the distribution of galaxies in redshift space, using redshift rather than true distance as the 
radial coordinate. In redshift space, peculiar velocities displace galaxies along the line of sight 
and introduce a preferred direction in the pattern of galaxy clustering. One can exploit this clus- 
tering anisotropy to measure properties of the galaxy velocity field and constrain the value of O, 
the cosmological density parameter (Sargent & Turner 1977). This paper describes some practical 
methods for carrying out such measurements, and it explores their limitations. 

In the linear regime of gravitational instability models, the anisotropy of clustering takes a very 
simple form in Fourier space, as shown by Kaiser (1987; see also McGill 1990). The strength of an 
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individual plane wave is amplified by a factor that depends on the angle between the wavevector 
and the line of sight, 



where and denote the Fourier amplitudes in real and redshift space respectively, and //ki 
is the cosine of the angle between the wavevector, k, and the line of sight, 1.^ The degree of 
amplification is controlled by /3 = where /(O) ~ O"'^ is the logarithmic derivative of the 

fiuctuation growth rate (see Peebles 1980, §14), and b is the "bias factor", an assumed constant 
of proportionality between galaxy and mass fiuctuations. In linear theory, (3 fixes the ratio of the 
peculiar velocity field to the gravitational acceleration field inferred from the galaxy distribution. 
Upon taking the ensemble average of the modulus of (1-1), one sees that the power spectrum in 
redshift space is anisotropic and given by 



where the real space power spectrum, P^(k), is assumed to be an isotropic function of k. Kaiser 
remarked that the angular dependence of the redshift space power spectrum could be used to 
measure the quantity /3, and hence to constrain O if the bias parameter could be determined 
independently. 

This paper extends Kaiser's analysis in several ways. First, we show that the expression (1.2) 
has a simple expansion in harmonic moments, and that the ratio of the quadrupole and monopole 
moments of P^(k, ji) provides a simple estimator for [i in the linear regime. Second, we examine the 
effects of non-linear evolution on P^{k,ii) using cosmological iV-body simulations. Non-linearities 
can be important out to quite large scales, so they complicate the task of inferring [i from the power 
spectrum, but they also split the degeneracy between O and h that exists in linear theory. Finally, 
we describe a practical method of measuring P^{k,ii) from galaxy redshift surveys. Equation (1.2) 
assumes that the observer is at a great distance from the volume occupied by the survey, so that 
all the galaxies in the survey can be treated as though they lie along the same line of sight. We 
develop a way to handle this constraint, given a redshift survey of (fairly) arbitrary geometry, and 
we show how to correct the estimated power spectrum for the systematic effect of shot noise and 
for the effect of selecting galaxies within a finite window. We test our method by applying it to 
mock redshift surveys drawn from iV-body simulations, and we present preliminary results from 
an analysis of the 1.2 Jy IRAS galaxy redshift survey (Fisher 1992; Fisher et al. 1993a; hereafter 



Many previous studies of clustering anisotropy have employed the redshift space correlation 
function, measured as a function of separations parallel and perpendicular to the line of sight. On 
small scales, the random velocities of galaxies in virialized clusters cause the well known "finger-of- 
God" effect, which stretches clusters along the line of sight in redshift space. On large scales, infiow 
compresses overdense regions along the line of sight, and outfiow expands underdense regions. The 
"finger-of-God" effect stretches contours of the correlation function along the line of sight at small 
separations, while streaming motions compress the contours at large separations. Davis & Peebles 
(1983) and Bean et al. (1983) used the small scale distortion to estimate the pairwise velocity 
dispersion of galaxies. The predicted form of the distortion on large scales was first calculated by 
Lilje & Efstathiou (1989), and later by McGill (1990), and Hamilton (1992). Hamilton (1993a) has 
estimated [i from the large scale anisotropy of the correlation function, using the IRAS 2 Jy survey 
of Strauss et al. (1992). Fisher et al. (1993b) have measured the redshift space correlation function 
of the IRAS 1.2 Jy survey and examined both small scale and large scale anisotropies. 

The correlation function and the power spectrum form a Fourier conjugate pair, so complete 
knowledge of one is equivalent to complete knowledge of the other. However, estimators of the 

^ Throughout this paper we adopt the notation /Uxixa the cosine of the angle between xi and 
X2, and we use the superscripts R and S to denote quantities measured in real and redshift space 
respectively. 
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two quantities may behave quite differently wlien applied to finite and noisy data sets. There are 
several advantages to characterizing large scale clustering in terms of the power spectrum instead of 
the correlation function. First, uncertainty in the mean galaxy density affects only the k = mode 
of the power spectrum, while standard techniques for measuring the correlation function suffer 
random and systematic errors at all scales if the mean density is not known a priori (see Hamilton 
[1993b] for partial solutions to these problems). Second, the error properties of the power spectrum 
are cleaner than those of the correlation function, which greatly simplifies the quantitative testing 
of models. In particular, if primordial fiuctuations are Gaussian, then Fourier modes of the density 
field are statistically independent in the linear regime, though one must take account of correlations 
introduced by estimating the power spectrum from a finite data sample. Third, theoretical studies 
suggest that the transition between Fourier modes that are evolving according to linear theory 
and shorter wavelength modes that have gone non-linear is fairly sharp (e.g. Little, Weinberg & 
Park 1991). Fourier transforming the spectrum mixes different physical scales, complicating the 
transition between the linear and non-linear regimes and hence the comparison between theory 
and observation. Finally, the power spectrum is more immediately related to the predictions of 
theoretical models, and on large scales the mapping from real space to redshift space is most 
naturally described in the Fourier regime. 

A number of authors have estimated (k) directly from galaxy redshift catalogues (Baumgart 
& Fry 1991; Gramman & Einasto 1991; Peacock 1991; Peacock & Nicholson 1991; Vogeley et al. 
1992; Park, Gott & da Costa 1992; FDSYH; Feldman, Kaiser & Peacock 1993). The analysis in 
this paper overlaps most with the studies of the IRAS redshift surveys by FDSYH and Feldman et 
al. (1993), and with the iV-body study of Gramman, Cen & Bahcall (1993). However, those papers 
examine only the spherically averaged power spectrum — the monopole moment of P^(k, ji) — and 
they concentrate on the shape of the spectrum as a function of k. We focus instead on determining 
O and h from the angular anisotropy of P^(k, ji). We also undertake a more systematic exploration 
of non-linear effects. 

Section 2 discusses theoretical aspects of the redshift space power spectrum. In §2.1 we derive 
a linear theory estimator for [i based on harmonic moments of the power spectrum. In §2.2 we 
explore non-linear effects using cosmological iV-body simulations. In §2.3 we relate our analysis 
to the correlation function analysis of Hamilton (1992). Section 3 presents a practical method for 
measuring the redshift space power spectrum and estimating [i from galaxy redshift surveys. In §4 
we apply this method to the 1.2 Jy IRAS redshift survey. Section 5 summarizes our conclusions. 
The two appendices contain mathematical details of the calculations in §2.3 and §3. Throughout 
the paper, we express distances in h~^Mpc, where h = i7o/(100 km s~^ Mpc~^). 

2 Theory 

2.1 A LINEAR THEORY ESTIMATOR FOR /3 

It is convenient and useful to characterize the redshift space power spectrum in terms of multi- 
pole moments. These moments define the decomposition of P^{k,ii) into Legendre polynomials, 
hereafter denoted Viiji), 



oo 




(2.1) 



1=0 



The multipole moments, Pf(k), can be computed by the inversion formula 



+1 




(2.2) 



-1 



For reference, the first three even Legendre polynomials are Vo(i-i) = 1, V2(l-i) = (3//^ — l)/2, and 
Viiji) = (35ii^ - 30//2 + 3)/8. 
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In linear theory, P^{k,ii) has the simple angular dependence of equation (1-2), which is char- 
acterized completely by its monopole, quadrupole, and hexadecapole moments. Equation (1.2) 
contains no odd powers of ji and no powers higher than so all odd multipoles and all moments 
with / > 4 must vanish according to equation (2.1). The non- vanishing moments can be obtained 
easily by direct integration: 



3^ 5 

4 „ 4 



Pi{k)= (i + -/3 + -/32)p«(fc) 



Pj;{k) = for / / 0, 2, or 4 . 

As we will see in §2.3, this behaviour is analogous to that of the redshift space correlation function 
described by Hamilton (1992). When clustering is non-linear, one can still perform a multipole 
decomposition of the power spectrum as in equation (2.1), but in this regime one expects that 
/ > 4 moments will contribute and that the low-order moments will depart from the values in 
equation (2.3). In particular, the "finger-of-God" distortion reverses the sign of the quadrupole 
term on small scales because it weakens fluctuations along the line of sight instead of amplifying 
them. 

If we have measured P^{k,ii) on large scales and decomposed it into multipole moments, then 
we can estimate [i from the ratio of the quadrupole term to the monopole term. The linear theory 
prediction, which follows directly from (2.3), is 

This function is well behaved over the interesting range of /3, growing almost linearly from G ~ 0.13 
at /3 = 0.1 to G ~ 1.02 at /3 = 1, so the inversion from G{p) to [i is stable and straightforward. 
In principle one can also obtain [i from the ratio Pl{k)l P^ik), but in practice we find that the 
hexadecapole is too sensitive to noise and non-linearities in P^{k,ii) to yield useful estimates. 

One could devise other schemes for extracting [i from measurements of P^(k, //), but there are 
two advantages to the estimator defined by equation (2.4). First, the monopole and quadrupole 
moments provide a "low-order" description of the anisotropy predicted by equation (1.2), so an 
estimator based on them should be relatively insensitive to noise. Second, the power spectrum 
estimator described in §3 is most naturally expressed in terms of multipoles, because the correction 
for a finite data sample is described most simply in terms of the correction to each multipole 
moment. It therefore makes sense to base the (3 estimator directly on these multipole moments. 

The techniques described by FDSYH and Feldman et al. (1993) yield estimates of the direction- 
averaged redshift space power spectrum, i.e. the monopole component of P^{k,ii). One can then 
derive [i by measuring the amplification of the redshift space power spectrum over the real space 
power spectrum, but only if the real space spectrum has itself been estimated by some other method 
(e.g. from an angular catalog). Our approach, on the other hand, exploits the angular dependence 
of the redshift space spectrum at fixed k, and the real space power spectrum cancels out of the 
estimator (2.4). 



2.2 NON-LINEAR EFFECTS 

In order to understand how non-linearities can affect P^{k,ii) and the [i estimator defined by 
equation (2.4), it is helpful to recap Kaiser's (1987) linear theory derivation of the power spectrum 
anisotropy, equation (1.2). For simplicity, we will assume a volume limited galaxy sample, so that 
terms involving derivatives of the selection function vanish. If the small scale velocity dispersion 
can be neglected on the scale of interest, then the galaxy density in redshift space is related to 
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the galaxy density in real space by the Jacobian of the transformation between the two coordinate 
systems, 



^ I U(r) - U(0) 
r 



-2 



\ , dU(r) 



dr 



2.5 



where U{y) is the radial component of the peculiar velocity. The first multiplying factor describes 
the effect of the peculiar velocity on the cross sectional area of an element of fixed angular size, 
while the second factor describes compression or expansion along the line of sight. The peculiar 
velocity associated with a plane wave of wavenumber k and amplitude A is f/ ~ A/A;, while 
dlJ/dr ~ klJ ~ A. We can therefore set the first factor to unity for a distant observer at r ^ 
obtaining 

p^(s) = p«(r)(l+^) ' • (2.6) 
If \dU/dr\ ^ 1, then an expansion of (2.6) yields 

= P"(r) (l - ^) ' (2-7) 
=^ S^s)^p^s)/p-l = S-ir)-^-S-ir)^ . (2.8) 

Now consider a single Fourier mode of the galaxy density field, with wavevector k and modulus 
Ak. The associated real-space density perturbation is 

^"(r) = AkCos(k • r + v') , (2.9) 

and in linear theory the corresponding velocity perturbation is 

v(r) = -/3 — Aksin(k-r + ¥)) . (2.10) 

In equation (2.10) we have assumed a linear bias between the galaxy and mass perturbations, 
Ak = bA'^"''^'^ . The radial component of the velocity perturbation is 

U(r) = v(r) • r = -^^Ak sin(k ■ r + Lp) , (2.11) 

=^ ^ = -Mki/3AkCos(k-r+¥') = -//2,;3^«(r) . (2.12) 

Combining this result with equation (2.8) yields 

^^(s) = S-(r)[l + + Mki/3^"(r)] (2.13) 
= S-(r)(l + plf^), if |^"(r)|<l. (2.14) 

The last result is simply equation (1.1), and equation (1.2) for the power spectrum follows upon 
taking the appropriate ensemble average. 

Four distinct guises of the linear approximation appear in this derivation, the assumptions 

that: 

(1) the small scale velocity dispersion can be neglected, 

(2) \dU/dr\ ^ 1, so that the Jacobian can be expanded as it is in equation (2.7), 

(3) the velocity and density perturbations satisfy the linear theory relation (2.10), 

(4) |^"(r)| ^ 1, so that equation (2.14) follows from equation (2.13). 

We have also made the distant observer approximation to obtain equation (2.6), but this is a 
geometrical approximation that is not specific to linear theory. Assumptions (1) — (4) all hold on 
scales that are truly in the linear regime. However, it is helpful to distinguish them because they 
break down at different rates as one moves to smaller scales. Furthermore, the relative severity of 
the breakdowns will vary from model to model, since the accuracy of assumption (1) depends on 
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the amplitude of small scale dispersions, (2) on the magnitude of large scale velocity gradients, (3) 
on the contrast of mass density fluctuations, and (4) on the contrast of galaxy density fluctuations. 
The multiplicity of non-linear effects makes their combined behaviour complicated, but it also offers 
the hope of achieving separate constraints on model parameters that would be degenerate in linear 
theory. For instance, in a model with very low O, one would expect violations of (f ) and (2) to be 
small relative to violations of (3) and (4). 

In what follows, we will refer to violations of assumption (I) as "dispersion non-linearity," 
violations of (2) as "gradient non-linearity," violations of (3) as "dynamical non-linearity," and 
violations of (4) as "contrast non-linearity." The boundary between dispersion and gradient non- 
linearity is not a clear one, since the very existence of dispersion in the velocity field contradicts 
the notion of a simple Jacobian transformation between real space and redshift space. Nonetheless, 
it is useful to maintain some conceptual distinction between a form of non-linearity that derives 
primarily from virialized objects and a form that refiects the amplitude of larger scale velocity 
gradients. Operationally, we can regard dispersion non-linearity as the component that is removed 
by smoothing the galaxy velocity distribution on small scales (a few Mpc). It is difficult to calculate 
the impact of these various non-linear effects analytically, even if they are treated in isolation. 
For now, we will limit ourselves to a preliminary investigation based on cosmological iV-body 
simulations. 

The iV-body simulations that we will analyze are a subset of those used by Little & Weinberg 
(1993) in their study of the void probability function. Each simulation has Gaussian initial condi- 
tions with the parametrized power spectrum of Efstathiou, Bond & White (1992). With Efstathiou 
et a/.'s scale parameter F set to 0.5, this power spectrum is a good fit to that of the standard 
cold dark matter model with 0=1 and h = 0.5. We adopt a lower value of the scale parameter, 
F = 0.25, to obtain a spectral shape closer to that implied by recent measurements of the power 
spectra and correlation functions of galaxies and galaxy clusters (Maddox et al. 1990; Dalton et 
al. 1992; Peacock & Nicholson I99I; Vogeley et al. 1992; FDSYH; Feldman et al. 1993). Each 
simulation occupies a periodic cube 300h~^Mpc on a side, and there are two independent realiza- 
tions of each model. We evolve the simulations into the non-linear regime using a staggered-mesh 
PM code written by Changbom Park (Park 1990), with 100"^ particles on a 200"^ density-potential 
mesh. In this paper we present results from models with 0=1 and O = 0.3, bias factors 5 = 1 
and 2, and cosmological constant A = 0. We have also examined models with O = O.I and with 
5 = 3. Initial power spectra are normalized so that the rms amplitude of mass fiuctuations in 
spheres of radius 8h~^Mpc would be cig = 1/5 if the initial conditions were extrapolated to 2; = 
by linear theory. For the unbiased (5 = I) models we select a "galaxy" population randomly from 
the iV-body particles. For the biased models, we select a biased subset of the particles according 
to the high peaks algorithm (Bardeen et al. 1986), so that for galaxies (t| = ha^ = I. More details 
of the simulations and the biasing procedure appear in Little & Weinberg (1993). From our point 
of view, the important feature of these simulations is their large physical volume, which allows 
us to measure the power spectrum with reasonable accuracy on large scales. However, the force 
resolution (I.5h~^Mpc grid scale) is rather low, which probably causes us to underestimate the 
magnitude of small scale velocity dispersions. The simulations are adequate for the exploratory 
purposes of this paper, but higher resolution simulations would be desirable for a more detailed 
examination of velocity dispersion effects. 

The light solid lines in Figure I show contours of the linear theory power spectrum (1.2) for 
(3 = 1/2 and P^(k) equal to the F = 0.25 power spectrum of our initial conditions. The axes fcn and 
kj_ indicate wavenumbers parallel and perpendicular to the line of sight, so //ki = ^ii/(^^ + k^Y^"^ ■ 
Note that the outer contours correspond to smaller spatial scales (higher k), and that the contour 
levels increase inwards. In linear theory, coherent fiows amplify fiuctuations with wavevectors close 
to the line of sight, so contours of the power spectrum are elongated in the fcn direction. If we 
imagine revolving these contours around the fcn axis (since there is a single line of sight direction 
but a full 27r range of tangential directions), we can describe this elongation as a prolate distortion 
of the spherically symmetric, real space power spectrum. 
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Figure 1. Contours of the redshift space power spectrum P {k, fi) in the plane defined by and kj_, the 

2 2 1 /2 

components of the wavevector parallel and perpendicular to the line of sight — = ^||/(^|| + ^j^) • Note 
that the innermost contours correspond to the largest physical scales. Contour levels increase from (k, /t) = 2 
at the outside to P^ (k, fi) = 32 at the inside. Light solid contours show the linear theory power spectrum, 
P^ (k, /t) = (1 + Pii^y^ P^(k), for /? = 1/2 and the T = 0.25 real space spectrum P^(k). Heavy solid contours 
show P^ (k, fi) measured from an iV-body simulation with f2 = 1 and 6 = 2; non-linear effects are substantial for 
all but the innermost contour. Dotted contours show the power spectrum of the same simulation after the particle 
velocities have been smoothed with a spherical filter of radius 5h ^ Mpc. Smoothing removes "fingers-of-God" and 
restores the contour excursions to large . 

The heavy contours in Figure 1 show the power spectrum measured from one of our 0=1, 
5 = 2 simulations. We shift the biased particles into "redshift" space, treating the the 2;— component 
of the peculiar velocity as the radial velocity, and we compute the power spectrum by Fast Fourier 
Transform (FFT). By identifying the 2;— axis with the line of sight, we automatically satisfy the 
distant observer constraint, and we ensure that the radial velocity field is periodic in the simulation 
volume, so that the boundary conditions assumed by the FFT are appropriate. Of course we 
cannot measure power spectra in the real universe this way; we leave that problem to §3. At 
the largest scale (innermost contour), the simulation power spectrum lies fairly close to the linear 
theory prediction, but at smaller scales there are substantial departures. In particular, the outer 
contours of the simulation power spectrum are oblate rather than prolate. 

The oblate shape of these outer contours is a signature of dispersion non-linearity, which washes 
out small scale fiuctuations along the line of sight. The reduction in small scale power compresses 
contours along the fcn axis. Oblate outer contours of P^(A;j^, fcn) are the Fourier analogue of prolate 
inner contours of the correlation function ^^(rp^ir), which are also caused by "fingers-of-God" 
(see Fisher et al. 1993b). The dotted contours of Figure 1 show the power spectrum of the same 
simulation, except that we have smoothed the velocity field on small scales by replacing the velocity 
of each particle with the mean velocity of all particles within a 5h~^Mpc sphere centred upon it. 
This smoothing procedure removes the velocity dispersion of virialized clusters and groups, but it 
leaves the large scale fiow intact. Smoothing alters the shapes of the outer contours dramatically, 
restoring their extensions to large k^^. Comparison of dotted and solid curves shows that dispersion 
non-linearity significantly distorts the shapes of all but the innermost contour. The two sets of 
contours always agree when k^Jk^ because peculiar velocities do not alter the power spectrum 
at //ki = 0. 

By squashing contours towards small k^^ (small //ki), dispersion non-linearity suppresses the 
quadrupole moment of P^{k,ii) relative to its monopole moment. It therefore depresses the value 
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Figure 2. Effect of dispersion non-iinearity on ttie iinear ttieory /? estimator (2.4). Ttie ratio of ttie estimated 
/? to ttie true value is plotted as a function of wavelength A = 2-k /k. a) A model with f2 = 1 and 6 = 2. Solid 
curves show results for two simulations without velocity smoothing; dispersion non-linearity depresses the estimated 
/? below the true value out to large scales. Dotted curves show results after smoothing the particle velocities over 
5h ^Mpc; the estimated /? rises rapidly to the true value. 6) Same as a), but for simulations with f2 = 0.3 and 
6 = 2. Smoothing has less effect in this model because velocity dispersions are smaller. 



of [i inferred by our linear theory estimator (2.4). Figure 2a illustrates this effect for the O = f , 
5 = 2 model. Solid lines show the ratio of the estimated [i to the true value, /3true = 1/b = 0.5, 
as a function of wavelength A = 2'iT/k. The two lines represent our two independent simulations; 
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Figure 3. Effect of gradient non-iinearity on ttie /? estimator (2.4). Soiid and dotted fines sfiow resufts for unbiased 
(b = 1) modeis witti f2 = f and f2 = 0.3, respectiveiy. Particie vefocities tiave been smootfied over 5fi ^Mpc in eacti 
simuiation, and vaiues pfotted are tfie average resufts from our two independent runs of eacfi modef. Tfie difference 
between tfie sofid and dotted curves reffects tfie fiigfier ampfitude of tfie farge scafe vefocity fiefd in tfie f2 = f 
modef. Tfie dasfied fine, nearfy coincident witfi tfie dotted fine, is estimated from tfie f2 = f modef after scafing aff 
of its particfe vefocities (and tfie vafue of /?true) by a factor of O.s"'*', to matcfi its vefocity ampfitude to tfiat of tfie 
n = 0.3 modef. 

the close match between them indicates that sample-to-sample fluctuations in P^(A;,//) are small 
for wavelengths A ^ 60h~^Mpc, in volumes the size of our 300h~^Mpc simulation cubes. The 
dotted lines in this flgure show the estimated values of /3//3true for the simulations after smoothing 
the velocity held. These curves rise much more quickly, leveling out at /3est//3true ~ 0.9 by A = 
40h~^Mpc. Comparison of the solid and dotted lines shows that dispersion non-linearity has a 
signiflcant impact in this model even at A = 60h-^Mpc. 

Figure 2b makes the same comparison for an O = 0.3, 5=2 model. Once again the curves 
for the velocity-smoothed simulations rise more rapidly to their asymptotic value, but the effect 
of velocity dispersions is much smaller in this low-0 model, as expected. Observational estimates 
of the galaxy pairwise velocity dispersion yield Gy ~ 350kms~^ at separations ~ lh~^Mpc (Davis 
& Peebles 1983; Fisher et al. 1993b; other estimates range from ~ 250kms~^ by Bean et al. 1983 
to ~ 450kms~^ by de Lapparent, Geller & Huchra 1988). The pairwise dispersion of the 0=1, 
5 = 2 simulations is 600kms~^, and that of the O = 0.3, 5=2 simulations is 250kms~^, so the 
impact of dispersion non-linearity in the real universe is probably smaller than that indicated by 
Figure 2a but larger than that of Figure 2b. 

The solid line in Figure 3 shows the [i inferred by our quadrupole estimator for the unbiased 
(5 = 1), O = 1 model, averaged over the two simulations. We have smoothed the velocity held 
over 5h~^Mpc as before, to remove the effects of dispersion non-linearity. Nonetheless, the /3est 
vs. A curve rises slowly, and it levels out at only /3est ~ ^-Tfitme- The dotted line shows the 
corresponding estimate for the unbiased, O = 0.3 model. Here the /3-estimator performs much 
better, rising to a value /3est ~ 0.85/3ti.ue by A = 20h~^Mpc, then leveling off. We have removed 
the small scale velocity dispersion by smoothing, and the large scale spatial clustering in the two 
models is virtually identical because they have the same initial power spectrum. The difference in 
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Figure 4. Effect of dynamicaf non-finearity on tfie /? estimator (2.4). Ttie soiid fine represents ttie modei witti 
n = 1 and 6 = 2 (/? = 0.5). Ttie dotted fine represents ttie modet witti = 0.3 and b = 1 (jS = 0.49). Botti 
curves are averaged over two simutations, and particte vetocities tiave been smoottied over 5ti ^Mpc. Ttie dastied 
tine stiows ttie dotted curve muttiptied by a factor of 1/0.85, so ttiat it tevets out at /?est//9true ~ 1- Atttiougti true 
l3 vatues are nearty identicat in ttie two modets, mass ftuctuations are non-tinear out to targer scates in ttie unbiased, 
tow-fi modet, so ttie /Sgst i^) curve tevets out at targer wavetengttis. Dynamicat non-tinearity offers a potentiat toot 
for obtaining separate constraints on f2 and b, given a sufficientty targe data set. 

results must therefore reflect the difference in the amplitude of the large scale velocity held, which 
influences /3est through gradient non-linearity. To conflrm this point, we scale down the individual 
particle velocities of the 0=1 model by a factor O.S"'^, then smooth the velocity held, measure 
the power spectrum, and estimate (3 as before. We obtain the dashed line of Figure 3, which is 
virtually indistinguishable from the dotted line that represents the O = 0.3 model. The detailed 
match of individual features in these curves arises because we use the same initial random phases 
for the = 1 and O = 0.3 simulations. 

Increasing (3 amplifles both dispersion and gradient non-linearity. We have already remarked 
that there is no clean separation between these two forms, but since we have smoothed away the 
small scale dispersions for both of the models shown in Figure 3, it seems appropriate to ascribe the 
remaining difference between them to gradient non-linearity. In the unbiased, = 1 model, this 
effect signiflcantly depresses the (3 inferred by our quadrupole estimator out to quite large scales. 
Results for an unbiased, O = 0.1 model are similar to those for the O = 0.3 model in Figure 3, 
which suggests that the residual impact of gradient non-linearity is small once O ^ 0.3 (/3 ^ 0.5), 
at least with this degree of velocity smoothing. 

In the linear regime, distortions of the redshift space power spectrum depend on a speciflc 
combination of physical parameters, /3 = O^'^/^. Clearly we would like to know the values of O and 
b individually. In particular, given the evidence that (3 < 1 from small scale clustering (e.g. cluster 
mass-to-light ratios or the cosmic virial theorem), we would like to be able to distinguish a low 
density universe in which galaxies trace mass from a critical density universe with biased galaxy 
formation. Can we exploit the scaling of non-linear effects to make such a distinction? 

To make the question speciflc, suppose that we analyze a large redshift sample and flnd that 
our quadrupole technique yields a stable estimate of /3 ~ 0.5 on large scales. We wish to distinguish 
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between two alternative hypotheses: an O = 0.3 model with 5=1 (and /3 = O.S"'^ = 0.49), and an 
= 1 model with 5 = 2. The amplitude of large scale flows in these two models is the same, and 
small scale velocity dispersions will be similar, so gradient and dispersion non-linearities cannot 
help us. The observed power spectrum directly determines the amplitude of galaxy fluctuations, 
and hence the importance of contrast non-linearity. However, the amplitude of mass fluctuations 
is different in the two models, with the unbiased, O = 0.3 model having larger fractional mass 
fluctuations at a given physical scale. One might therefore expect that dynamical non-linearity 
would appear in the low-0 model on a scale where the 0=1 model still satisfles equation (2.10), 
the linear theory relation between velocity and density. 

Figure 4 illustrates just this behaviour. The solid line shows the ratio /3est//3true for our = 1, 
5 = 2 simulations, with the velocity held smoothed over 5h~^Mpc. It rises to /3est ~ Arue by 
A ~ 12h~^Mpc, then levels off. The dotted line represents the O = 0.3 model, and it does not level 
off until A ~ 20h~^Mpc. The "turnover wavelength" of the unbiased model is nearly double that of 
the biased model, reflecting the importance of dynamical non-linearity at larger scales. In the region 
where /3est is stable, A ~ 20 — 40h~^Mpc, the value for the O = 0.3 model is only /3est ~ 0.85/3ti.ue- 
The difference from the = 1 model, which has /3est ~ Arue in this region, presumably reflects 
the continuing influence of dynamical non-linearity, since other non-linear effects should be similar 
in the two models. With a sufficiently large observational sample, we would expect to recover 
Pest ~ /3true ^t larger wavelengths. However, presented with the dotted curve and no a priori 
knowledge of O, we might mistake the stable value found between 20 and 40h~^Mpc for /Jtrue, 
and rescale the height of the /3est//3true vs. A curve accordingly. The result would be the dashed 
curve of Figure 8, which is simply the dotted curve multiplied by 1/0.85. Although this curve lies 
closer to the solid curve that represents the = 1 results, one can still discern the larger scale of 
non-linearity in the unbiased model. 

The moral of this section is that non-linear effects are a nuisance, but a potentially useful 
nuisance. Non-linearity can affect the anisotropy of the redshift space power spectrum out to 
quite large scales — wavelengths of 60h~^Mpc or more — and the variety of effects prevents a 
simple analytic treatment. We have catalogued these effects and obtained some sense of their 
numerical importance, but we have not yet devised practical means for handling or exploiting 
them. Dispersion non-linearity seems likely to pose the thorniest problems. We have dealt with 
it here by smoothing the velocities in real space, but that approach is not possible with a redshift 
catalog. One strategy is to smooth the data in redshift space (equivalent to multiplying in the 
Fourier domain), perhaps using an anisotropic fllter that takes advantage of the known anisotropy 
of "flngers-of-God." Another option is to estimate the pairwise velocity distribution on small 
scales from the redshift space correlation function, then use this distribution function to correct 
the power spectrum. Of course one can also compare the uncorrected data directly to theoretical 
models that include small scale dispersion, but here one confronts the problem that the pairwise 
velocity dispersion is a difficult quantity to compute reliably from theory (witness the continuing 
controversy over "velocity bias"). 

The dependence of dynamical non-linearity on mass fluctuations rather than galaxy fluctua- 
tions offers the tantalizing prospect of constraining O and the bias factor separately. Unfortunately, 
dispersion, gradient, and contrast non-linearities can mask these dynamical effects, making them 
difficult to exploit. Further analytic and numerical investigations should help this effort; a dis- 
tinctive signature that can isolate dynamical non-linearity from other effects would be especially 
useful. The distinction between the unbiased, low-0 model and the biased, = 1 model in Figure 
4 is subtle, difficult to measure even in our 300h~^Mpc cubes of "perfect" data. From the current 
generation of galaxy redshift surveys, an estimate of (3 is probably the most we should expect. How- 
ever, huge redshift samples like the million-galaxy Sloan Digital Sky Survey (Knapp & Gunn 1993) 
will allow high precision measurements of the galaxy power spectrum, so separate determinations 
of O and 5 should be possible. 
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2.3 COMPARISON TO CORRELATION FUNCTION ANALYSIS 

It is informative to compare the method of estimating (3 described in §2.1 to Hamilton's (1992; 
1993a) method, which is based on the redshift space correlation function ^^(s,!^). Hamilton de- 
composes the redshift space correlation function into harmonic moments ■^f (s), defined by 

e(s,ii) = Y,^i('^)nj^)- (2.15) 

He then estimates the rational function G(P) [equation (2.4)] from a rather elaborate combination 
of the quadrupole and monopole terms, 

^^^^ = -^a-^) + (Js^) j:m:^'y'd.s' ' ^'-''^ 

and derives (3 accordingly. 

In Appendix B, we show that the harmonic components of ^^(k,!^) can be expressed elegantly 
in terms of the harmonic components of P^(k, ji). We give expressions for these components and for 
some linear combinations of (s), such as the denominator of (2.16), which are useful in comparing 
our formalism to that of Hamilton (1992). Here we make use of these identities to re-express (2.16) 
in terms of the Pj^ (k). Substituting the identities (B.5) and (B.8) into (2.16), we obtain a simple 
and informative expression for Hamilton's estimator of G(P), 

. 1^ dkk'^P^{k)i2{ks) 

^^^^ = j^dkk^P^{k)Uks) ' ^^-^^^ 

where j2{kr) is the usual spherical Bessel function. In the linear regime the ratio P2(k)/PQ(k) is 
constant and equal to G(P), and equation (2.16) therefore yields an unbiased estimate of (3. The 
effects of non-linearity enter at large k, where P2ik) is depressed or even changes sign and Pq (k) 
is boosted. These effects depress the value of G(P) given by (2.16) at small s. Non-linear small 
scale power can affect G(P) out to quite large separations because of the oscillating weight function 
i2{ks). If the transition between linear and non-linear behaviour is sharp in the power spectrum, 
it will be fuzzy in the correlation function. 

The power spectrum has relatively simple error properties on large scales, which give it some 
practical advantage over the correlation function. However, the crucial question is which quantity 
can be estimated more accurately on scales in the linear regime, and the answer to this question 
may be survey- and model-dependent. 



3 Measuring P^{k,ii) from a redshift survey 



Kaiser's formula for the redshift space power spectrum, equation (1.2), is valid only if is 
computed from data that subtend a small solid angle on the sky. Without this constraint the 
quantity //ki is not even clearly defined, since there is no single line of sight direction 1. Faced with 
the task of measuring P^{k,ii) and estimating [3 from a galaxy redshift survey, our approach is to 
first extract a subset of the redshift data such that all of the selected galaxies lie approximately 
along the same line of sight. Specifically, we select galaxies within a spherically symmetric window 
of radius Rg-ph centred at a distance L from the observer sufficient that 9 ~ L/Rsph ^ 1- Only 
galaxies that fall inside the window are used in estimating the power spectrum. However, we 
can repeatedly reposition the window within the volume defined by the geometry of the redshift 
survey and the requirement that L ^ Rsph, then average the results to obtain a final, combined 
estimate that utilizes most of the galaxies in the survey. A roving window of this form has practical 
advantages, since one can apply the same type of analysis to cone surveys or to "all sky surveys" 
with incomplete coverage near the galactic plane. In practice, the limited depth of existing redshift 
surveys forces one to compromise between using distant, poorly sampled data, which satisfy ^ ^ 1 
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but yield a noisy Fourier transform, or carrying out the transform with nearby, well sampled 
structure, for which the small-angle approximation is less accurate. At the end of this section we 
will calibrate the effect of a finite opening angle on the power spectrum estimate empirically, using 
a catalogue drawn from an iV-body simulation. 

Once the data are windowed to ensure an acceptably small opening angle, the estimates of 
Fourier modes at different k are no longer independent. Selecting galaxies within a window is 
equivalent to multiplying the galaxy density field by a window function (which is zero outside the 
windowed region). In Fourier space, this multiplication becomes a convolution with the Fourier 
transform of the window function, which introduces correlations between neighbouring modes over 
a range Ak ~ 1/Rsph- In any real application, one therefore expects the measured P^{k,ii) to 
represent the true redshift space power spectrum convolved with a function related to the Fourier 
transform of the window used to define the sample. This convolution tends to make the measured 
redshift space power spectrum more isotropic, and therefore, if uncorrected, it causes a system- 
atic underestimate of [i. Fortunately, this effect can be calculated analytically (see below and 
Appendix A), and corrections can be applied to the estimated power spectrum. 

Suppose that we wish to measure the power spectrum from a redshift survey characterized by 
mean galaxy density (ra) and radial selection function (^(r) — i.e., in the absence of inhomogeneities 
the mean galaxy density at distance r would be {n)(j){r). A straightforward estimator for the Fourier 
mode amplitudes is the sum 

N 



{n)V^ i=i 



where N is the total number of galaxies in the survey and Vyj = j d^s w{s) is the weighted volume 
of the sample. The summation in equation (3.1) is carried out in redshift space, and windowing 
of the data is introduced explictly through the window function w{s). Again, we have in mind a 
window function that is non-zero over a region (e.g., a sphere) whose origin is displaced far enough 
from the central observer that it subtends a small angle. The simplest such function is a "top-hat" 
defined by centre and radius Rg-ph'- w{s) = — sd/Rsph), where is a unit step function. The 
sharp edges of a top hat create nasty sidelobes in its Fourier transform, so in our analysis we adopt 
a "bowler-hat" window, defined as the convolution of a top hat with a Gaussian: 



VJ{ s 



(2nr'^'Rf I d\'^\s'-s,\/R,p,)e-\^-^'\y'^l (3.2) 



Gaussian convolution apodizes the edges of the window and reduces the sidelobes in its Fourier 
transform. Results in this paper were computed using Rg/Rsph = 0.1. 

We show in Appendix A that the ensemble average of the squared-modulus of 5k is 

{S^S^) = P''''ik,i,u) + V^\wik)f + I d^s^-^^ , (3.3) 

where w(k) = V,^~^ J d'^s w(s)e~^^^'^ is the Fourier transform of window function, normalized to 
unity at k = 0. The second term on the right hand side of (3.3) refiects the finite width of the 
window function in Fourier space; in the limit of a large volume, this term would contribute only 
at k = 0. The last term represents the discreteness or "shot noise" contribution to the Fourier 
transform. The term of interest is P^^ (k, i^i^i), which is the convolution of the true redshift space 
power spectrum P^(A;,//ki) with the square of the window transform: 

P^-(A;,Mki) = ^ I d'k'P'{k',fi^.,)\w{k-k')f . (3.4) 



Multipole decomposition allows a simple treatment of the convolved redshift space power 
spectrum. If we define Pj^^(k) to be the multipole moments of the convolved spectrum P^^ {k, ji), 
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Convolution Kernel 




k'/k 

Figure 5. Tlie convolution kernels ^i(k, k') that appear in the convolution of the multipole moments, equa- 
tion (3.5), computed for our bowler-hat window function. The solid curve refers to the monopole moment (l = 0), 
while the dashed and dotted curves refer to the quadrupole (l = 2) and hexadecapole (l = 4) moments, respectively. 



analogous to the decomposition of {k^ji) in equation (2.1), we can derive a relation between the 
true and convolved multipole moments (cf.. Appendix A): 

oo 

Pf'/k) = j dk'e'Pf{k')Ti{k,k') , 



^ pn.mik) . 

Equation (A. 15) of Appendix A defines the convolution kernel !Fi{k,k'). The "transfer" functions 
{k) describe the effect of the convolution on the multipole moments. In general, they depend on 
the wavenumber k, on the form of the power spectrum Pf(k), and on the window geometry. How- 
ever, for sensible choices of the window function w(s), the convolution kernel !Fi{k,k') is strongly 
peaked around k = k' , and in this case the integral that defines (k) depends only very weakly 
on k and the power spectrum. 

For a given value of k, we scale the radius of the window function Rg-ph to give w{k) = 0. The 
Fourier transform of a Gaussian is positive definite, so for the bowler-hat window the first zero 
occurs at the first zero of the top hat, which is Rg-ph = 4:A934/k. This scaling has the advantage 
of making the second term in equation (3.3) identically zero. The only remaining term to be 
subtracted from {S\^S^) is independent of k, so the shape of the estimated power spectrum (in 
a linear plot) is insensitive to uncertainties in the mean density (ra) (FDSYH). This choice also 
simplifies the form of the transfer functions. The window transform w(|k — k'|) is a function of 
the combination |k — k'|i?sp/i, so with Rgph scaled to k it becomes a function of |k — k'|/A; = 
[1 -|- {k' jky -\- 2//kk'^'/^]^^^- Substituting into equation (A. 15) and integrating over //kk', we find 



(3.5) 
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Figure 6. Recovery of the redshift space power spectrum by our "practical" algorithm. Format is similar to that 
of Figure 1. Heavy contours, repeated from that figure, show the redshift space power spectrum measured directly 
by FFT from an f2 = 1, 6 = 2, iV-body simulation. Dotted contours show the power spectrum recovered by the 
algorithm described in §3, which can also be applied to galaxy redshift surveys. These contours are reconstructed 
from just the monopole, quadrupole, and hexadecapole components. They agree well with contours of the FFT 
measurement at large scales. At small scales (outer contours), non-linearity introduces some power into the higher 
harmonics of P^'ik, fi), so the recovery is imperfect. 

that the convolution kernels !Fi{k,k') depend only on the ratio k' jk. If P^(k) is scale-free, then 
the functions {k) are independent of k. 

Figure 5 shows the convolution kernels !Fi{k,k') for our bowler-hat window, as functions of 
k' jk. The solid curve refers to the monopole moment (/ = 0), while the dashed and dotted curves 
refer to the quadrupole (/ = 2) and hexadecapole (/ = 4) moments respectively. The multipole 
moments of the redshift space power spectrum measured in a finite window are the true multipole 
moments convolved according to equation (3.5) with the curves shown in Figure 5. 

Table 1 

Spectral Index T2 /Tq T^/Tq 



n = 


-2 


0.680 


0.333 


n = 


-1 


0.706 


0.384 


n = 





0.716 


0.422 


n = 


+ 1 


0.711 


0.445 



Table 1 gives the values of the transfer coefficients T^"^ for our bowler-hat window, in the case 
where the true spectrum multipoles are Pfik) oc k^. From Table 1, it is clear that the effect 
of convolution on the multipole moments of P^{k,ii) is insensitive to the shape of the real space 
power spectrum. This is good news, since it means that one can correct the observed multipole 
moments of P^^ {k^ji) to recover the true redshift space power spectrum P^ {k^ji) without a priori 
knowledge of the power spectrum. However, Table 1 also shows that the convolution is a large 
effect for our choice of window function, so it cannot be neglected. For example, the ratio of the 
measured quadrupole to the measured monopole is suppressed by 30% from its true value. 

We can test our "practical" method of estimating P^ {k^ji) by applying it to one of the iV-body 
simulations described in §2.2. We first select a value of A; and define a bowler-hat window function 
with Rs-ph = 4:A934/k and Rg = O.lRsph- We view the simulation along a randomly chosen line of 
sight, 1, and locate the bowler-hat window randomly within the simulation volume. For now, we 
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Figure 7. The value of /? inferred from tire linear theory estimator (2.4). The solid line shows estimates based 
on the FFT measurement of (k, fi) averaged over the two f2 = 1, 6 = 2, iV-body simulations. Circles show the 
result of applying our "practical" algorithm to the same simulations. Results of the two methods agree well, though 
the estimated /9's lie below the true value /? = 1/2 because of the continuing effects of dispersion non-linearity on 
these scales (see Figure 2). 

define radial velocities as if we were viewing the data from an infinite distance, i.e. in the limit of 
a small opening angle. We carry out the Fourier sum 5k, defined by equation (3.1), for a set of 
50 randomly oriented k-vectors. We average the squared-moduli S\^S^ in bins of A//ki = 0.1 and 
repeat the whole process iV^, times, choosing iV^, > L\^oxl R^gph to ensure that we adequately sample 
all the structure in the simulation. Subtracting the shot noise term of equation (3.3) leaves us with 
an estimate of P^^ {k, ji). We fit its angular dependence using a basis of Legendre polynomials 
Viiji) and hence determine the coefficients Pf^'(k) in its multipole expansion [equation (A. 8)]. 
To reconstruct the unconvolved P^{k,ii) we need the multipole moments Pf(k), which we obtain 
from the Pf^(k) using equation (3.5) and the transfer coefficients T^^ {k) computed for our window 
function and the F = 0.25 power spectrum. 

The heavy contours in Figure 6 show the redshift space power spectrum of the O = 1,5 = 2 
iV-body simulation, measured directly by FFT, just as in Figure 1. The dotted contours show 
the power spectrum measured as above, reconstructed from just the monopole, quadrupole, and 
hexadecapole terms of P^{k,ii). The two measurements of the power spectrum agree well at large 
scales. At the small scales represented by the outermost contours, non-linearity has introduced 
power into the higher harmonics of P^(A;,//), so our three-term reconstruction does not yield a 
complete description. By including more terms of the harmonic expansion, we could presumably 
recover P^{k,ii) more accurately in this regime. 

We do not need a complete reconstruction of P^{k,ii) to apply our [i estimator (2.4), just 
the monopole and quadrupole moments. The solid line of Figure 7 shows the value of [i estimated 
from the direct FFT measurement of P^{k,ii), as in Figure 2a. Open circles show the estimates 
of [i based on the monopole and quadrupole moments recovered by our "practical" method. The 
two estimates agree well at all wavelengths. However, the estimated /3's lie below the true value, 
especially at smaller wavelengths, because of the non-linear effects discussed in §2.2, with dispersion 
non-linearity being the most serious of these. To obtain accurate estimates of [i from a real redshift 
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Figure 8. The ratio Rq of tlie quadrupole measured wlien tire windowed data subtend a finite opening angle 8 to 
tfie value obtained for a vanishingly small opening angle. This estimate was made for a wavelength A = 80Mpc/h 
in an f2 = f, 6 = 2 iV-body simulation. However, the effect is purely geometrical, and therefore independent of 
wavelength or the intrinsic value of /?. Corrections for the finite opening angle must be applied when 8 exceeds ~ 30 
degrees. 

survey, one must either correct for non-linearities or measure the power spectrum on still larger 
scales. 



In a real redshift survey we cannot place our window at infinite distance, so the data sample 
must subtend a finite opening angle. When the opening angle is large, the power spectrum estimator 
(3.3) averages over a range of line- of- sight directions //ki at fixed k, which depresses the amplitude 
of the measured quadrupole moment and thereby reduces the estimated value of (3. We can use the 
iV-body simulations to examine the magnitude of this effect. To do so, we measure P^{k,ii) with 
our bowler-hat window function, but we set the observer at a distance L = Rg-phl tan(^/2) from the 
window centre. Figure 8 shows the ratio Rq of the quadrupole measured with opening angle 9 to 
the asymptotic value obtained with a vanishingly small opening angle, for < ^ < 90°. We see that 
Rq > 0.98 for 9 ^ 30°, i.e. the measured quadrupole is diluted by less than 2% and the resulting 
estimate of (3 is nearly unbiased. The dilution effect increases rapidly at larger angles. We have used 
a wavelength A = 80h~^Mpc to compute Figure 8, but the dilution is a purely geometrical effect 
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reflecting the spread in //ki values within the window, so the quantity Rq should be independent 
of wavelength and of the intrinsic value of [i. When analyzing a real redshift survey we can either 
ensure that the opening angle is less than 30°, or, if the limited survey depth requires that we use 
a larger opening angle, we can make a systematic correction to the measured quadrupole based on 
the ratio Rq. 

Our algorithm for measuring the redshift space power spectrum and estimating [i can be 
summarized as follows: 

(1) For each wavenumber k, pick a random line of sight 1 and place a bowler-hat window at a 
distance L > R^ph/ ^^T^iOmax/'^-), so that it subtends an opening angle 9 < 9max- The window 
radius Rgph = 4.4934/A; is chosen to make the Fourier transform = 0. 

(2) For many random choices of the wavevector direction k, carry out the Fourier sums (3.1) over 
galaxies in the window to estimate the mode amplitudes 5k. Compute {S\^S^), the squared- 
moduli averaged in bins of //kb subtract the shot noise term of equation (3.3) to obtain 
an estimate of P^'^(A;, //ki), the convolved redshift space power spectrum. 

(3) Increase the signal-to-noise ratio by repeating this procedure for many random choices of the 
window centre, within the volume deflned by the survey boundaries and the opening angle 
constraint, and averaging the results. 

(4) Fit the multipole moments Pf^'(k) of the measured P^'^(A;, //ki). Correct these using the 
transfer functions {k), deflned by equation (3.5), to obtain the moments Pi {k) of the true 
redshift space power spectrum. If the maximum opening angle O^ax exceeds 30°, divide the 
quadrupole moment by the geometrical correction factor Rq plotted in Figure 8, and apply 
similar corrections for any higher moments of interest. 

At this point, one can reconstruct P^{k,ii) from its estimated multipole moments Pf(k), and one 
can use the ratio P2(k)/PQ(k) to estimate (3 via equation (2.4). The (3 estimator is valid only in 
the linear regime, but the technique for measuring P^{k,ii) is general and applies on all scales. 

4 Preliminary Application to the IRAS 1.2 Jy Survey 



We now present a preliminary application of our method to the 1.2 Jy IRAS survey studied 
by FDSYH. The sample we analyse is complete to a flux limit of 1.2 Jy at 60 //m and contains 
5304 galaxies with galactic latitude |5| > 5° . A detailed description of the sample and its selection 
criteria can be found in Strauss et al. (1990) and Fisher (1992). 

We use the method described in §3 to measure the angular dependence of the redshift space 
power spectrum, P^(A;,//), at a variety of wavelengths. For each wavelength A = 27r/A;, we use 
bowler-hat data windows [equation (3.2)], with Rg-ph = 4:A934/k and Rg = OARsph- We centre 
each window at a random angular position on a shell of radius L = iZ^p/i/ tan(55°/2), avoiding any 
overlap with the sample boundary at |5| = 5°. The longer the wavelength the more distant the shell, 
and therefore the larger the contribution of shot noise to the measured power. We adopt the rather 
large opening angle of ^ = 55° in order to keep the shot noise tolerably small. The large opening 
angle slightly reduces the amplitude of the measured quadrupole component of the power spectrum, 
as described in §3. We correct this geometrical effect using the appropriate value of Rq read from 
Figure 8. In order to sample the data on each shell completely, we take 500 sub-samples deflned 
by different, random angular positions of the bowler-hat window. In each sample we estimate the 
power for 50 randomly oriented k-vectors, using the selection function (^(s) adopted by FDSYH. 
We average the resulting estimates as a function of //kb i^i bins of width A//ki = 0.1. 

At wavelengths A = 25,30,35 and 40h~^Mpc, the estimated power spectrum shows a clear 
enhancement for wavevectors close to the line of sight. The estimate of P^{k,ii) at A = 50h~^Mpc 
is considerably noisier than at the shorter wavelengths. This limitation is inevitable with the 
present data set, since at 50h-iMpc the shot noise term that is subtracted from the initial power 
estimate is as large as the remaining power. Figure 9 plots the values of [i deduced from the ratios 
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Figure 9. Values of /? estimated from tlie redsliift space power spectrum of tire 1.2 Jy IRAS redshift survey, using 
the linear theory estimator (2.4), as a function of wavelength A = 2Tr/k. Results for 30h ^Mpc< A < 45h ^ Mpc 
are consistent with /? = 0.35 it 0.05. Non-linear effects remain important on these scales, so the true value of /? is 
probably higher. Errors due to shot noise become large beyond A ~ 45h ^ Mpc. 

of the quadrupole to monopole components at each of these wavelengths, after applying the transfer 
function correction appropriate for spectral index n = —1 (from Table 1) and the 4% geometrical 
correction for the 55° opening angle (from Figure 8). In the range 30h~^Mpc< A < 45h~^Mpc, 
the estimates are consistent with (3 = 0.35 ± 0.05. However, our experience with the iV-body 
simulations shows that the effect of velocity dispersions in groups and clusters is not negligible on 
these scales, and that this value is likely to be an underestimate. 

This analysis of the 1.2 Jy IRAS survey is preliminary in several senses. Most obvious is 
the fact that Figure 9 contains no error bars. An estimate of the errors under the assumption of 
Gaussian fluctuations could be carried out using the formalism of Feldman et al. (1993), but a more 
thorough error analysis would require the use of mock IRAS catalogues constructed from iV-body 
simulations, along the lines of FDSYH. We plan to perform such an analysis in a future, detailed 
study of the 1.2 Jy survey. Coupled to the absence of error bars is our decision to use, at each 
wavenumber k, only that portion of the survey that lies on a shell defined by the opening angle 
9 = 55°. Utilizing the whole survey would require us to combine estimates from different shells with 
different shot noise, and taking a proper weighted average of the results would require knowledge of 
the errors on individual estimates. Because the shot noise increases rapidly with distance, it is not 
clear that we would gain much by including more distant shells in our power spectrum estimates, 
but a modest improvement in signal-to-noise should be possible. Finally we have not made any 
correction for non-linear effects, in particular for the effects of small scale velocity dispersions. We 
plan to use the simulated IRAS catalogues mentioned above to investigate possible methods of 
supressing dispersion non-linearity, such as compressing clusters or smoothing the data in redshift 
space. By combining results from different wavelengths, utilizing the survey volume in a more 
optimal manner, and correcting for dispersion non-linearity, we should be able to improve our 
estimate of (3 significantly. 
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5 Conclusions 

By shifting galaxies along the line of sight, peculiar motions make the redshift space power 
spectrum of galaxy clustering, P^(A;,//), anisotropic. In linear perturbation theory, the degree of 
anisotropy depends only on the parameter (3 ~ O^'^/^, and the form is characterized entirely by 
its monopole, quadrupole, and hexadecapole moments. One can estimate the value of (3 from the 
ratio of the quadrupole and monopole moments of P^{k,ii), for any wavenumber k in the linear 
regime. 

We have described and tested a practical method of measuring P^{k,ii) from galaxy redshift 
surveys, which uses an offset data window to select galaxies that lie close to the same line of sight. 
This method applies in both the linear and non-linear regimes. Key features of the method are a 
correction for discreteness contributions to the power spectrum and a correction for the systematic 
effect of selecting galaxies within a finite region. There are many advantages to analyzing redshift 
space distortions using the power spectrum. The prediction of linear theory takes a very simple form 
in Fourier space. Estimates of the power at different wavenumbers are statistically independent, 
provided they are separated by more than the width of the window in Fourier space, so estimates 
of [i at different wavelengths can be combined in a straightforward way to yield a more accurate 
overall estimate. Also, estimates of the power spectrum are insensitive to errors in the mean density 
of a survey, unlike estimates of the correlation function on large scales. 

We have used a series of iV-body simulations to explore the effects of non-linear gravitational 
evolution on the power spectrum anisotropy. Non-linearities take a variety of forms, and they can 
infiuence the power spectrum out to wavelengths of 50h~^Mpc or more. These effects complicate the 
task of inferring [i from P^{k,ii) — typically, they cause our quadrupole method to underestimate 
the true value of [i. However, the transition between the linear and non-linear regimes depends 
separately on O and 5, so with sufficiently good data it should be possible to constrain both 
parameters individually. 

The present study leaves at least three challenging tasks for future work. The first is to devise 
a practical method for removing or minimizing the effects of small scale velocity dispersions, since 
these seem to be the most pernicious and least useful form of non-linearity. The second is to 
calculate the uncertainties and correlations of power spectrum estimates and to use the results of 
this calculation to develop an optimal strategy for analyzing fiux-limited data sets and combining 
estimates of [i from different scales. The third is to obtain a more thorough understanding of 
"dynamical" non-linearity — the breakdown of the linear theory relation between velocity and 
density perturbations — since this form of non-linearity is the most useful for separating O and h. 

Our preliminary analysis of the 1.2 Jy IRAS redshift survey reveals a clear signal of anisotropy 
in the power spectrum on large scales. The quadrupole-to-monopole ratio yields (3est = 0.38, 0.32, 
and 0.42 at wavelengths A = 30, 35, and 40h~^Mpc, respectively. The residual effects of small 
scale velocity dispersions remain important at these wavelengths, so we expect these values to 
underestimate the true value of (3. The limited depth of the 1.2 Jy survey prevents us from making 
reliable estimates at significantly larger scales, but suitable corrections for non-linearity should 
allow us to make the estimate from these wavelengths more accurate. 

Prospects for applying this technique to other data sets are encouraging. The sparsely sampled 
QDOT survey of IRAS galaxies goes to a fainter fiux limit than the 1.2 Jy survey, and it will soon 
be extended to a complete survey of IRAS galaxies brighter than 0.6 Jy. On the optical side, the 
extension of the CfA survey is nearly complete, extension of the Southern Sky Redshift Survey is 
underway, and the multi-fibre survey of Kirshner, Oemler, Schechter & Schectman is progressing 
rapidly. All of these data sets should provide useful, and largely independent, estimates of (3. The 
million-galaxy redshift survey of the Sloane Digital Sky Survey should yield a precise measurement 
of (3 and separate constraints on O and b. 
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Appendix A. The Convolved Redshift Space Power Spectrum 

Suppose that we have a redshift survey characterized by a mean galaxy density (ra) and radial 
selection function <?!>(r), defined so that the number of galaxies per unit volume at distance r would 
be {n)(j){r) in the absence of inhomogeneities. As an estimator for the Fourier amplitudes, consider 
the quantity 

(n)V,S S 

where N is the total number of galaxies in the survey and Vyj = j d?sw{s). The window function 
w{s) defines the limits of the data sample used to compute 5k, since data lying in the region where 
w{s) = do not contribute to the Fourier sum. However, w{s) need not be constant over its non- 
zero region, and it need not be centred on the observer. In this paper, for example, we adopt the 
bowler-hat window described in §3, and we centre it far enough from the observer that it subtends 
an opening angle 9 < 55°. 

Following Peebles (1980, §41.1), we rewrite 5k as a sum over infinitesimal cells with volumes 
d'^Si and occupation numbers A^: 

= E . (A.2I 

{n)V.l,t!:, 

The ensemble average of the squared-modulus of 5k is 

The expected occupation numbers of cells located at positions and Sj are 

For the i = j term, we have used the fact that Ni must equal or 1 in infinitesimal cells, implying 
Nf = Ni and (Nf) = (Ni). In real space {S^(xi)S^*(xj)) depends only on the distance x = |xi — Xj| 
and is given by the usual two-point correlation function, ^^(x). In redshift space the correlation 
function, like the power spectrum, is anisotropic (Kaiser 1987; Lilje & Efstathiou 1989; McQill 
1990; Hamilton 1992), and it depends both on the length of the separation vector, s = |si — Sj\, 
and on the vector's orientation with respect to the line of sight. The correlation function ^^(s^i^si) 
and the power spectrum P^(A;', //k'l) form a Fourier transform conjugate pair, with 



(^^(s,)^^*(s,))^e'(5,Msi) 



1^ I d'k'P'(k',,ik'i)e-''''-% (A.5) 



where s = Si — Sj. Substituting equation (A.5) into equation (A. 4) allows equation (A.3) to be 
written, after some manipulation, as 

('5k5i:) = -J% [ d^k'P^k',i,^,,)\wik-k')f+V^\wik)f + j^ [d^s^-^^, (A.6) 

(2TTy J {n)Vuj J ns) 
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where w(k) = V^^ J vj(s)e~^^^'^d'^s is the Fourier transform of the window function, normalized to 
unity at k = 0. 

The first term in equation (A. 6) is the true redshift space power spectrum convolved with the 
squared Fourier transform of the data window. The second term arises from the finite width of the 
window function in Fourier space; in the limit of a large volume, this term would contribute only 
at k = 0. The last term is the discreteness or "shot noise" contribution to the Fourier transform. 

When estimating the power spectrum for a given value of A;, we choose the radius of the 
data window to yield w(k) = (following FDSYH). We then compute the squared-modulus of the 
Fourier sum defined by equation (A.l) and subtract the shot noise (final) term of equation (A. 6). 
The line of sight is the only preferred direction, so we can average results over different directions 
of k that have the same //kb we can average over different positions of the data window within 
the redshift survey. We are left with an estimate of the convolved redshift space power spectrum. 



(A.7) 



The effect of the convolution on P^{k,ii) can be conveniently expressed in terms of multipole 
moments. We proceed by first decomposing the true and convolved power spectra and the window 
transform factor \w{k — k')p into their multipole moments, 

oo oo 

P'{k,^^) = Y^p^kmn) , P'^k,^^) = Y^pr{k)Vi{n) , 

oo 

and \w(k - k')|' = Y. Mk, k')Vi{i^kk,) . (A.8) 

1=0 

The expansion coefficients can be obtained from the inversions 

+1 +1 

pf{k) ^ ^ I di,p^k,i,)Vi{i,) , prik) ^ ^ I di,p'-^\k,i,)Vi{i,) , 
-1 -1 

+1 

wi{k, k') = I di^^^, \w (k - k')|' P,(Mkk') . (A.9) 

-1 

In equations (A.8) and (A.9), Viiji) denotes the Legendre polynomial of order /. 

With the expansions given in (A.8), the convolved power spectrum appearing in (A.7) can be 
written as 

oo 

P"^(A;,M = EE / f^fc'fc"^f(fc')*K^,^')/f^f^k'7',(Mk'i)7'KMkk') . (A.IO) 

I' I Q 

We can carry out the sum over /' by applying the following theorem for Legendre polynomials, 

j dn^, P,(Mk'sjPKMk'sJ = ^^Viil^s.sMi' , (A.ll) 
where Sfi, is the Kronecker delta symbol. Theorem (A.ll) reduces equation (A.IO) to 



P'"(A;,Mki) = E 



oo 

'^''^ ' -J lJ2 vs nJ\„r. 



21+1 (27r)3 





dk' k'^Pf(k')wi(k,k') 



Viifiki) . (A.12) 
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Recalling the definition of wi{k,k'), this becomes 

+1 



I 



{2 



^ j dk'k'^Pfik') j dfi^^, \w{k-k')fVi{fikk, 



kl 



(A.13) 



By comparing (A.13) with the expansion for P^^{k,ii) given in (A. 8), we see that the terms in the 
square brackets are just the multipole moments of the convolved spectrum, Pj^^(k), and that these 
are related to the unconvolved coefficients by the following transfer function, 



pf(k) 



1 



Pf(k)(2 



oo +1 

^ I dk'k'^Pfik') I di^^^, \w{k-k')fVi{i^kk, 



(A.14) 



Pfik) 



dk' k'^Pf{k')Ti{k,k') , 







where the convolution kernel appearing in the last line is 

+1 

Hk,k' 

-1 



(A.15) 



Note that J-i(k, k') depends only on the window geometry, not on the power spectrum. Furthermore, 
for a sensible choice of window !Fi{k,k') is sharply peaked around k = k' , so the power spectrum 
factors inside and outside the integral of equation (A.14) nearly cancel, making the transfer function 
insensitive to the shape of Pf {k) (cf. Table 1). In the limit that the window function w{s) encom- 
passes a large volume, the effect of the convolution become negligible, Ti{k,k') k~^S^^\k — k'), 
and we recover {k) = 1. 



Appendix B. The Redshift Space Correlation Function 

In this Appendix we derive identities relating the harmonic components of the redshift space 
correlation function to those of the redshift space power spectrum. We use these identities to show 
that in the linear regime our formalism reproduces the expressions derived by Hamilton (1992). 

The redshift space correlation function and power spectrum form a Fourier conjugate pair, 

e'(5,Msi) = ^ j d'kP'{k,fi^,)e-^>^^>^-^ . (B.l) 

The vector 1 denotes the line of sight direction, and //ks denotes the cosine of the angle between 
the vectors s and k. We now replace P^{k,ii) by its expansion in terms of Legendre polynomials, 

P'{k.^i,.^) = Y,Pnk)Vl{^i^) , (B.2) 

I 

and we substitute for e~''^*^ks using the Rayleigh expansion of a plane wave, 

e-'''^"^' = E(-0'(2/ + l)ji(ks)Vi(Hk,) , (B.3) 

where ji(ks) is the spherical Bessel function. We obtain 

oo 

ef(5,Msi) = EE(-0''(2/' + l)7^ / dkk'Pf(k)ji,(ks) i dnMm)Vi,(jik,) . (B.4) 
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Applying theorem (A. 11) yields the relation 



oo 



Clis) = I dkePf{k)ji{ks) , (B.5) 



where ^f(s) is the order-/ harmonic coefficient of the redshift space correlation function, i.e. 

r(5,Msi) = Eef(5)nMsi). (b.g) 

The odd-/ components of P^{k,ii) and ^^{s^fi) vanish by symmetry, so the factor i' is simply ±1 
and determines the relative sign of harmonic distortions in the correlation function and the power 
spectrum. The expansion defined by (B.6) and (B.5) is an identity, valid in both the linear and 
non-linear regime. In the case of linear theory, one can substitute for Pi {k) using the relations 
(2.3). The resulting harmonic expansion of ^^(s,!^) is identical to the expression given by Fisher 
et al. (1993c). 

In order to compare the Fourier analysis to the correlation function analysis presented by 
Hamilton (1992), the following identities are useful: 

oo 

^ j dkk'P^(k)jo(ks) = ^^(s) , (B.7) 



oo s 

^ I dkePi{k)3,{ks) = -eo'(^)+ ^ / ds's'^^^is') , (B.8) 



oo s s 

^ j dkePi{k)3,{ks) = a{s) + \^ j ds's''a{s')-\^ j ds's'^ais'). (B.9) 



The identity (B.7) is simply (B.5) for the case of / = 0, and the other two identities can be derived 
from (B.7) using the recurrence relations and integral properties of spherical Bessel functions. If 
one then substitutes the linear theory relations (2.3), these relations give equations (6), (7) and (8) 
of Hamilton (1992). 
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